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We analyze the perturbative and parametric stability of the QCD predictions for the Callan-Gross 
ratio R(x,Q 2 ) = Fl/Ft in heavy-quark leptoproduction. We consider the radiative corrections 
to the dominant photon-gluon fusion mechanism. In various kinematic regions, the following con- 
tributions are investigated: exact NLO results at low and moderate Q 2 < m 2 , asymptotic NLO 
predictions at high Q 2 3> m 2 , and both NLO and NNLO soft-gluon (or threshold) corrections at 
large Bjorken x. Our analysis shows that large radiative corrections to the structure functions 
Ft{x, Q 2 ) and Fl(x, Q 2 ) cancel each other in their ratio R(x, Q 2 ) with good accuracy. As a result, 
the NLO contributions to the Callan-Gross ratio are less than 10% in a wide region of the variables 
x and Q 2 . We provide compact LO predictions for R(x,Q 2 ) in the case of low i<1. A simple 
formula connecting the high-energy behavior of the Callan-Gross ratio and \ow-x asymptotics of the 
gluon density is derived. It is shown that the obtained hadron- level predictions for R(x — » 0, Q 2 ) 
are stable under the DGLAP evolution of the gluon distribution function. Our analytic results 
simplify the extraction of the structure functions F2(x,Q 2 ) and F2(x,Q 2 ) from measurements of 
the corresponding reduced cross sections, in particular at DESY HERA. 
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I. INTRODUCTION 

In the framework of perturbative quantum chromodynamics (QCD), the basic spin-averaged characteristics of 
heavy- flavor hadro- [l|, photo- Q and electro-production 0] are known exactly up to the next-to-leading order 
(NLO). Although these explicit results are widely used at present for a phenomenological description of available 
data (for reviews, see Refs. 0, Q), the key question remains open: How to test the applicability of QCD at fixed 
order to heavy-quark production? The basic theoretical problem is that the NLO corrections are sizeable; they 
increase the leading-order (LO) predictions for both charm and bottom production cross sections by approximately 
a factor of two. Moreover, soft-gluon resummation of the threshold Sudakov logarithms indicates that higher-order 
contributions can also be substantial. (For reviews, see Refs. @, Q-) O n the other hand, perturbative instability 
leads to a high sensitivity of the theoretical calculations to standard uncertainties in the input QCD parameters. For 
this reason, it is difficult to compare pQCD results for spin-averaged cross sections with experimental data directly, 
without additional assumptions. The total uncertainties associated with the unknown values of the heavy-quark mass, 
m, the factorization and renormalization scales, /if and the asymptotic scale parameter Aqcd and the parton 
distribution functions (PDFs) are so large that one can only estimate the order of magnitude of the pQCD predictions 
for charm production cross sections in the entire energy range from the fixed-target experiments [8j to the RHIC 
collider 

Since these production cross sections have such slowly converging perturbative expansions, it is of special interest 



'Electronic address: 
t Electronic address: 



nikiv@mail.yerphi.am 
kniehl@desy.de 



2 



to study those observables that are well-defined in pQCD. A nontrivial example of such an observable was proposed 
in Refs. [9L fiol fllL \l2\ fl3l| . where the azimuthal cos(2c/?) asymmetry in heavy-quark photo- and leptoproduction was 
analyzed. 1 ' 2 In particular, the Born-level results were considered [9( and the NLO soft-gluon corrections to the basic 
mechanism, photon-gluon fusion, were calculated (Io| . It was shown that, contrary to the production cross sections, 
the azimuthal asymmetry in heavy-flavor photo- and leptoproduction is quantitatively well defined in pQCD: the 
contribution of the dominant photon-gluon fusion mechanism to the asymmetry is stable, both parametrically and 
perturbatively. Therefore, measurements of this asymmetry should provide a useful test of pQCD. As was shown in 
Ref. [Til ], the azimuthal asymmetry in open charm photoproduction could be measured with an accuracy of about 
ten percent in the approved E160/E161 experiments at SLAC [15| using the inclusive spectra of secondary (decay) 
leptons. 

In Ref. [l3|], the photon- (heavy-)quark scattering contribution to (^-dependent lepton-hadron deep-inelastic scat- 
tering (DIS) was investigated. It turned out that, contrary to the basic photon-gluon fusion component, the quark- 
scattering mechanism is practically cos(2y)-independent. This is due to the fact that the quark-scattering contribution 
to the cos(2ip) asymmetry is, for kinematic reasons, absent at LO and is negligibly small at NLO, of the order of 1%. 
This indicates that the azimuthal distributions in charm leptoproduction could be a good probe of the charm PDF 
in the proton. 

In the present paper, we continue the studies of perturbatively stable observables by considering the photon-gluon 
fusion mechanism in heavy-quark leptoproduction, 

£(l) + N(p) -> 1(1 -q) + Q( PQ ) + X[Q](p x ). (1) 

In the case of unpolarized initial states and neglecting the contribution of Z-boson exchange, the cross section of 
reaction Q can be written as 



d 2 aiN Air a 



2 

cm 



dxdQ 2 Q 4 
2ira 2 



xQ 4 



{ [1 + (1 - y) 2 } F T (x, Q 2 ) + 2 (1 - y) F L (x, Q 2 )} 

{[! + (!- V?} F 2 (x, Q 2 ) - 2xy 2 F L (x, Q 2 )} , (2) 



where a cm is Sommerfeld's fine-structure constant, F2(x,Q 2 ) = 2x(Ft + Fl) and the kinematic variables are defined 
by 



S={i+p)\ Q 2 = -q 2 , 



Q 2 



2p-q' 
2 



v = — t , Q 2 = xyS, t=%. en 

p ■ t Tfl 

In this paper, we investigate radiative corrections to the Callan-Gross ratio in heavy-quark leptoproduction, defined 
as 

R ^ = WY (4) 

First, we consider the exact NLO corrections to the quantity R(x, Q 2 ) at low and moderate Q 2 < m 2 using explicit 
results [H, [T^] . Then, we analyze the high-Q 2 regime with the help of the asymptotic NLO predictions for the structure 
functions Ft(x,Q 2 ) and Fl(x,Q 2 ) presented in Refs. [13, [T^]. Finally, the soft-gluon (or threshold) contributions 
are investigated in the large- a; region in the framework of the formalism developed in Ref. Q. To next-to- leading 
logarithmic (NLL) accuracy, we calculate the NLO and NNLO soft-gluon corrections to both structure functions. Our 
main results can be formulated as follows: 

• Exact NLO corrections to the ratio R(x, Q 2 ) do not exceed 10% in the energy range x > 10~ 4 at low and 
moderate Q 2 < m 2 . 



1 Well-known examples include the shapes of differential cross sections of heavy-flavor production, which are sufficiently stable under 
radiative corrections. 

2 Note also the recent paper [13 ], where the perturbative stability of the QCD predictions for the charge asymmetry in top-quark 
hadroproduction has been observed. 
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• At high Q 2 ^> to 2 , the asymptotic NLO corrections to R(x, Q 2 ) are less than 10% for 10 4 < x < 10 1 . 

• At the NLL level, the NLO and NNLO soft-gluon predictions for R(x, Q 2 ) affect the LO results by less than a 
few percent at low and moderate Q 2 and x > 10~ 2 . 

• In all the cases mentioned above, the NLO predictions for R(x, Q 2 ) are sufficiently insensitive, to within ten 
percent, to standard uncertainties in the QCD input parameters p,p , (1r and Aqcd, and in the gluon PDF 
g(x,fi F ). 

We conclude that, in contrast to the production cross sections, the Callan-Gross ratio in heavy-quark leptoproduction 
is an observable quantitatively well defined in pQCD. Pcrturbative stability of the photon-gluon fusion results for 
R(x,Q 2 ) is mainly due to the cancellation of large radiative corrections to the structure functions Ft(x,Q 2 ) and 
Fl{x, Q 2 ) in their ratio, especially in the large- x region. Measurements of the quantity R(x, Q 2 ) in charm and bottom 
leptoproduction should provide a good test of the conventional parton model based on pQCD. 

Concerning the experimental aspects, perturbative stability of the QCD predictions for R(x,Q 2 ) observed in our 
studies is very useful for the extraction of the structure functions F%(x, Q 2 ) and F| (x, Q 2 ) from the data. Usually, it 
is the so-called "reduced cross section", a(x,Q 2 ), that can directly be measured in DIS experiments: 

~i rf\ 1 xQ 4 d 2 a [N 2 2xy 2 2 

°( x >Q) = TTTi 32 0—2 TZJ7v> = F 2 (x,Q)-—— -^F L (x,Q ) (5) 



1 + (1 - y) 2 2ira 2 m dxdQ 2 <v>-*; l + {l-y) 2 

= F 2 {x 7 Q 2 ) 



"~ R2(x,Q 2 ) 



i + (i - y) 



2 J 



(6) 



where 



P i pp. \ 9 F L (x,Q 2 ) R(x,Q 2 ) 

R ^ Q ] = 2x J^m = l + R(x,Q 2 Y (7) 

In earlier HERA analyses of charm and bottom electroproduction [l9l |. the corresponding longitudinal structure 
functions were taken to be zero for simplicity. In this case, a(x,Q 2 ) = F 2 (x,Q 2 ). In recent papers [13, HH> the 
structure function F 2 (x, Q 2 ) is evaluated from the reduced cross section (jS} where the longitudinal structure function 
Fl{x,Q 2 ) is estimated from the NLO QCD expectations. Instead of this rather cumbersome procedure, we propose 
to use the expression ^ with the quantity R 2 (x, Q 2 ) calculated in LO approximation. This simplifies the extraction 
of F 2 (x, Q 2 ) from measurements of a(x, Q 2 ) but does not affect the accuracy of the result in practice. 

Indeed, the LO corrections to the extracted function F 2 (x, Q 2 ) due to the non-zero value of R 2 (x, Q 2 ) cannot exceed 
30% because the ratio R 2 (x, Q 2 ) is itself less than 0.3 practically in the entire region of the variables x and Q 2 . For 
this reason, the NLO corrections to R 2 (x,Q 2 ), having a relative size of the order of 10%, cannot affect the value 
of F 2 (x, Q 2 ) by more than 3%. In reality, the effect of radiative corrections to R 2 (x, Q 2 ) on the extracted values of 
F 2 (x,Q 2 ) is less than 1% since y <C 1 in most of the experimentally accessible kinematic range. 

In the present paper, we derive compact hadron-level LO predictions for the ratio R 2 {x,Q 2 ) in the limit of low 
x — ► 0. Assuming the low- a; asymptotic behavior of the gluon PDF to be of the type g(x, Q 2 ) oc l/x 1+s , we provide 
analytic result for the ratio R 2 {x — ► 0, Q 2 ) = R^\Q 2 ) for arbitrary values of the parameter 6 in terms of the Gauss 

hypergeometric function. Furthermore, we consider compact formulae for R^\Q 2 ) in two particular cases: (5 = 1/2 
and (5 = 0. The simplest case, 5 — 0, which has already been studied recently in Ref. [22j| . leads to a non-singular 
behavior of the structure functions for x — > 0. The second choice, 6 = 1/2, historically originates from the BFKL 
resummation of the leading powers of ln(l/a;) 23]. 

In principle, the parameter S is a function of Q 2 and this dependence is calculated using the DGLAP evolution 
equations [24[. However, our analysis shows that hadron-level predictions for R 2 (x — > 0, Q 2 ) depend weakly on 
(5 practically in the entire region of Q 2 for (5 > 0.2. In particular, the relative difference between R^'^^Q 2 ) and 

R 2 3 \Q 2 ) is l ess than few percent at Q 2 m 2 . For this reason, our simple formula for R^\q 2 ) with 6 = 1/2 (i.e., 
without any evolution) describes with good accuracy the low-x predictions for R 2 (x, Q 2 ) of the CTEQ PDF versions 
[25l [26|. We see that the hadron-level predictions for R 2 (x — > 0, Q 2 ) are stable not only under the NLO corrections 
to the partonic cross sections, but also under the DGLAP evolution of the gluon PDF. 

Finally, we show that our compact LO formulae for R 2 S \Q 2 ) conveniently reproduce the HERA results for F% (x, Q 2 ) 
and F%(x, Q 2 ) obtained by HI Collaboration [13,[2l[ with the help of more cumbersome NLO estimations of Fl(x, Q 2 )- 

This paper is organized as follows. In Section HH we analyze the exact NLO results for the Callan-Gross ratio at 
low and moderate Q 2 < m 2 and the asymptotic NLO predictions at high Q 2 ^S> to 2 . The soft-gluon contributions to 
R(x,Q 2 ) are investigated in Section ITTT1 To NLL accuracy, we calculate the threshold NLO and NNLO corrections 
to both structure functions Ft(x,Q 2 ) and Fl(x,Q 2 ). The analytic LO results for the ratio R 2 (x,Q 2 ) at low x are 
discussed in Section HVl 
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FIG. 1: Feynman diagrams of photon-gluon fusion at LO. 

II. NLO PREDICTIONS FOR THE CALLAN-GROSS RATIO R(x, Q 2 ) 
A. Born-Level Cross Sections 

At LO, 0(a cm a s ), leptoproduction of heavy flavors proceeds through the photon-gluon fusion mechanism, 

l*(q)+g(k g )^Q( PQ ) + Q( PQ ). (8) 

The relevant Feynman diagrams are depicted in Fig. [1] The ^y*g cross sections, a^f 1 (z, A) (k — 2, T, L), have the form 

a£\z,\) = ^i^a B (z)|[(l-z) 2 + 2 2 + 4Az(l-3z)-8A 2 z 2 ]lni±|i-[l + 4z(l-z)(A-2)] / 3 z |, (9) 

4°)( 2 , A ) = H^M)fr fl ( z )J-2A«lni±^ + (1 -*)&}, (10) 
4 0) (z,A) = a ( °\z,\)-a { °\z,\), (11) 

with 



(2^) 2 e 2 a 



0111 



*B(Z) = — gT Z ' ( 12 ) 

where eQ is the electric charge of quark Q in units of the positron charge and a s {[i 2 R ) is the strong-coupling constant. 
In Eqs. ([§))- (fL2")l , we use the following definition of partonic kinematic variables: 



*=2?v A =<r ^-r-i-z (i3) 



Q' 2 x m 2 /„ 4Az 

!« 

The hadron-level cross sections, er^x, Q 2 ) (fc = 2, T, L), have the form 

l 

cr k (x,Q 2 )= J dzg(z,fi F )a k (^,X,fi F: fi R ^ : (14) 

x(l+4X) 

where g(z,HF) is the gluon PDF of the proton. The leptoproduction cross sections ak(x,Q 2 ) are related to the 
structure functions Fk(x,Q 2 ) as follows: 



F k (x,Q 2 ) = ® a k (x,Q 2 ) (k — T, L), (15) 



F 2 (x 1 Q 2 ) = — Q—a 2 (x,Q 2 ). (16) 

47T O cm 
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FIG. 2: Left panel: x dependence of the Callan-Gross ratio, R(x,Q 2 ) = Fl/Ft, in charm leptoproduction for £ = 0.44, 4.4 
and 44 at LO (solid lines) and NLO (dashed lines). Right panel: x dependence of the K factor for the transverse structure 
function, K(x, Q 2 ) = F^ LO /F^°, at the same values of £. 



B. Exact NLO Predictions at Low and Moderate Q 2 

At NLO, 0{a cni a 2 s ), the contribution of the photon-gluon component is usually presented in terms of the dimen- 
sionless coefficient functions c/™' l \z,\) (k = T,L), as 



, , e 2 Q a cm a s {n 2 j i 
a k {z,X,m ,n ) = 5 <^ c K k '(z,X) +47ra s (^ ) 



C --(z,A)+c--(z,A)ln^ 



+ 0(a 2 s ). (17) 



where we identify fi = /if = jiR. 

In this paper, we neglect the 7*<?(<?) fusion subprocesscs. This is justified as their contributions to heavy-quark 
leptoproduction vanish at LO and are small at NLO [3(. To be precise, the light-quark-initiated corrections to both 
Ft and Fu structure functions are negative and less than 10% in a wide kinematic range Our estimates show that 
these contributions cancel in the ratio R(x, Q 2 ) = Fl/Ft with an accuracy less than few percent. 

The coefficients (z, A) and {z, A) of the /Lt-dependent logarithms can be evaluated explicitly using renormal- 

ization group arguments The results of direct calculations of the coefficient functions c^'°'(z, A) (fc = T, L) are 

presented in Refs. [3l.[lq]. Using these NLO predictions, we compute the x dependence of the ratio R(x, Q 2 ) — Fl/Ft 
at several values of £ = 1/A = Q 2 /m 2 . 

The left panel of Fig. [2] shows the Callan-Gross ratio R(x, Q 2 ) as a function of x for £ = 0.44, 4.4 and 44. In our 
calculations, we use the CTEQ5M parametrization of the gluon PDF together with the values m c = 1.3 GeV and 
A3 = 373 MeV [26]. Unless otherwise stated, we use \x = ^/Am^~Q 2 throughout this paper. 

For comparison, the right panel of Fig. [5] shows the x dependence of the QCD correction factor for the transverse 
structure function, K(x,Q 2 ) = F^ LO /F^°. One can see that large radiative corrections to the structure functions 
Ft{x,Q 2 ) and Fl(x,Q 2 ), especially at non-small x, cancel each other in their ratio R(x,Q 2 ) = Fl/Ft with good 
accuracy. As a result, the NLO contributions to the ratio R(x, Q 2 ) are less than 10% for x > 10~ 4 at low and moderate 
Q 2 <m 2 c . 

Another remarkable property of the Callan-Gross ratio closely related to fast perturbative convergence is its para- 
metric stability. 3 Our analysis shows that the fixed-order predictions for the ratio R(x,Q 2 ) are less sensitive to stan- 
dard uncertainties in the QCD input parameters than the corresponding ones for the production cross sections. For 
instance, sufficiently above the production threshold, changes of \x in the range (1/2)^/4to 2 _j_ Q2 < ^ < 2y/4m 2 . + Q 2 
only lead to 10% variations of R(x, Q 2 ) at NLO. For comparison, at x = 0.1 and £ = 4.4, such changes of [i affect the 
NLO predictions for the quantities Ft(x,Q 2 ) and R(x,Q 2 ) in charm leptoproduction by more than 100% and less 
than 10%, respectively. 



Of course, parametric stability of the fixed-order results does not imply a fast convergence of the corresponding series. However, a fast 
convergent series must be parametrically stable. In particular, it must exhibit feeble fip and [ir dependences. 
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FIG. 3: Q 2 dependence of the asymptotic high-Q 2 (Q 2 3> m 2 ) predictions for the Callan-Gross ratio, R(x,Q 2 ) = Fl/Ft, in 
charm leptoproduction at x = 10 _1 , 10 -2 , 10 -3 and 10 -4 in LO (solid lines) and NLO (dashed lines). 



Keeping the value of the variable Q 2 fixed, we analyze the dependence of the pQCD predictions on the uncertainties 
in the heavy-quark mass. Sufficiently above the production threshold, i.e. in the plateau regions of the x distributions 
of R(x, Q 2 ) in Fig.0 changes of the charm-quark mass in the interval 1.3 GeV < m c < 1.7 GeV affect the Callan-Gross 
ratio by 2%-3% at Q 2 = 10 GeV 2 . The corresponding variations of the structure functions Ft(x, Q 2 ) and Fl(x, Q 2 ) 
are about 20%. We also verify that all the CTEQ versions [HI, H(| of the gluon PDF lead to NLO predictions for 
R(x, Q 2 ) that coincide with each other with an accuracy of about 5% practically at all x > 10~ 4 . 



C. Asymptotic NLO Results at High Q 2 » m 2 

The analytic form of the heavy-quark coefficient functions for lepton-hadron DIS in the kinematical regime Q 2 m 2 
is presented in Refs. (l7l . [l8j . The calculations were performed up to NLO in a s using operator product expansion 
techniques. 4 In the asymptotic regime £ — > oo, the production cross sections have the following decomposition in 
terms of the coefficient functions a!^ n ' m \z) (k — 2,L): 



- I 2 2\ 



l—l m-\-n<l 



o 



(18) 



It was found in Refs. [Tt], [3] that the hadron-level structure function F^ symp (x,Q 2 ) approaches, to within ten 

2 



percent, the corresponding exact value F| xact (x, Q 2 ) for £ > 10 and x < 10 1 both at LO and NLO. In the case 



For the longitudinal cross section Sx(z, Q ,m ,/x ), the asymptotic heavy-quark coefficient functions a L 
in a s O. 



I ,(n,m) 



(z) are known up to NNLO 
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FIG. 4: Q 2 -dependence of the asymptotic high-Q 2 (Q 2 S> m 2 ) predictions for the K factor, K(x,Q 2 
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F$ LU /F^ U , at 



of the longitudinal structure function F^ yinp (x, Q 2 ), the approach to F^ cact (x, Q 2 ) starts at much larger values of 
£ > 4 x 111-. 

Using the analytic NLO results for the coefficient functions presented in Ref. [l7|, we calculate the asymptotic 
high-Q 2 behavior of the ratio R(x,Q 2 ) — Fl/Ft at several values of x. Figure [3] shows R asymp (x,Q 2 ) in charm 
leptoproduction as a function of £ for x = 1CP 1 , 1CP 2 , 10~ 3 and 1CP 4 . In Fig. |4j we show the Q 2 dependence of 
the asymptotic predictions for the K factor K(x,Q 2 ) — F^ LO /F^° at the same values of x. One can see that the 
quantity K(x, Q 2 ) is practically independent of Q 2 at fixed values of x and tends to unity at low x. This implies that 
perturbative stability of the Callan-Gross ratio at low x is due to the smallness of the radiative corrections to both 
structure functions. At non-small x, the radiative corrections to Ft(x,Q 2 ) and Fl(x,Q 2 ) are large but cancel each 
other in their ratio R(x, Q 2 ) — Fl/Ft with good accuracy. 



III. SOFT-GLUON CORRECTIONS AT NLO AND NNLO 



In this Section, we consider the NLO and NNLO predictions for the Callan-Gross ratio due to the contribution of 
the photon-gluon fusion mechanism in the soft-gluon approximation and propose an improvement. For the reader's 
convenience, we collect the final results for the parton-level cross sections to NLL accuracy. More details may be 
found in Refs. @, LI3l • 

At NLO, photon-gluon fusion receives contributions from the virtual 0(ct em a 2 ) corrections to the Born process © 
and from real-gluon emission, 

7 *(g) + g(k g ) Q( PQ ) + Q(p Q ) + g(p g ). (19) 
The partonic invariants describing the single-particle inclusive (1PI) kinematics are 

s' = 2q ■ k g = s + Q 2 = (S", <i = (kg - p Q ) 2 - m 2 = (Ti, 

s 4 = s' + ti+Ui, ui = (q - pq) 2 - m 2 = Ui, (20) 

where £ is defined through k g = (p and S4 measures the inelasticity of the reaction (fl9|) . The corresponding 1PI 
hadron-level variables describing the reaction ([T]) are 

S' = 2q-p = S + Q 2 , Ti = (p - p Q ) 2 - to 2 , 

S 4 - S' + Ti + Ui, U x = (q-p Q f -m 2 . (21) 

The exact NLO calculations of unpolarized heavy-quark production in 73 [2|, 7*5 [3|, and gg [l| collisions show 
that, near the partonic threshold, a strong logarithmic enhancement of the cross sections takes place in the collinear, 
\p g ,T\ —> 0, and soft, \p g \ — * 0, limits. This threshold (or soft-gluon) enhancement is of universal nature in perturbation 
theory and originates from an incomplete cancellation of the soft and collinear singularities between the loop and the 
bremsstrahlung contributions. Large leading and next-to-leading threshold logarithms can be resummed to all orders 
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of the perturbative expansion using the appropriate evolution equations [29j . The analytic results for the resummed 
cross sections are ill-defined due to the Landau pole in the coupling constant a s . However, if one considers the 
obtained expressions as generating functionals and re-expands them at fixed order in a s , no divergences associated 
with the Landau pole are encountered. 

Soft-gluon resummation for the photon-gluon fusion was performed in Ref. 0] and confirmed in Refs. [ToL H3|. To 
NLL accuracy, the perturbative expansion for the partonic cross sections, d 2 ak{s' , t%, ui)/(dti dui) (k = T,L), can be 
written in factorized form as 



,2 d 2 <7 fc 



dtidui 

with the Born-level distributions Bf orn (s', t\, Ui) given by 



6(s' + h + Ui) + £ (^A " K (n) {s > i h)Ul) 



n=l 



i? orn (s', 



B 



Born// 



{- 


+ ir +4 (, 








f s m s' 







s'{m 2 -Q 2 /2) | Q 2 
Uui s' 



(22) 

(23) 
(24) 



Note that the functions K^ n \s' , ti,Ui) in Eq. (|22[) originate from the collinear and soft limits and are the same for 
both cross sections &t and &l- At NLO and NNLO, the soft-gluon corrections to NLL accuracy in the MS scheme 
read 



K^(s',ti,ui) = 2 



In (s 4 /m 2 ) 




' 1 " 






s 4 


+ 
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i 2 2C A 

where 62 = (HCU — 2n/) /12 is the first coefficient of the beta function, 

ding, (/j 2 ) ~ , a 

^ (Qs)= dinM 2 =-^ bk+ A~ 



k=l 



(26) 



(27) 



in Eqs. (|2"5|) and ([2"6"|) . CU = iV c , CV = (A*" 2 — 1)/(2N C ), rif is the number of active quark flavors, iV c is the number of 
quark colors, and Lp = (1 — 2r7i 2 /s){ln[(l — (3 Z )/(1 + {3 z )]+iir} with (3 Z = v/l — 4m 2 /s. The single-particle inclusive 
"plus" distributions are defined by 



In' (s 4 /m 2 ) 


= hm 


s 4 


+ 



In' (s 4 /m 2 ) 



0(s 4 - e) 



hi 



l+i 



S4, ■ I + 1 

For any sufficiently regular test function h(s 4 ), Eq. (|28p implies that 

In' (s 4 /to 2 ) 



r<K*4) 



(28) 



J dsi h(si) 




V (s 4 /m 2 )" 


s 4 


s 4 





ds 4 [/i(s 4 ) - ft(0)] 



,s 4 



i + 1 



/i(0) In 



i+i a i 



m - 



(29) 
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FIG. 5: Left panel: LO (solid lines), NLO (dashed lines) and NNLO (dotted lines) soft-gluon predictions for the x dependence 
of the Callan-Gross ratio, R(x,Q 2 ) = Fl/Ft, in charm leptoproduction at £ = 1, 2 and 5. Right panel: x dependence of the 
K factors K {1) {x,Q 2 ) = F^ LO /F^° (solid line) and K {2) (x,Q 2 ) = F^ NLO /F^ LO (dashed curve) for the transverse structure 
function at the same values of £. 



In Eqs. (P5|) and ([2rJ|) . we have also preserved the NLL terms for the scale-dependent logarithms. Note that 
Eqs. ([25] ) -ip5 ]) agree to NLL accuracy with the exact 0{a eva a 2 ) calculations of the photon-gluon cross sections or 
and <7l given in Ref. 

Numerical investigation of the results (p?5)) - (Prj)) was performed in Refs. @, It was shown that soft-gluon 

corrections reproduce satisfactorily the threshold behavior of the available exact results for the partonic cross section 
(72 = &t + cfl at £ < 1. Since the gluon PDF supports just the threshold region, the soft-gluon contribution dominates 
the hadron-level structure function at energies not so far from the production threshold. It was shown in Ref. Q 
that Eqs. (piS)) and ([^5]) render it possible to describe with good accuracy the exact NLO predictions [3| for the function 
F2(x,Q 2 ) at x > 1CP 3 and relatively low virtuality Q 2 ~ to 2 . 

In the present paper, we analyze separately the partonic cross sections ot and &l- It turns out that the quality 
of the adopted soft-gluon approximation is worse for than for &t- To clarify the situation, let us remember that 
the NLL approximation allows us to determine unambiguously only the singular S4 behavior of the cross sections 
defined by Eq. This implies that the S4 dependence of the Born-level distributions Bj,°£ n (s', t\, 1*1) \ u _ s _ s ,_ t 

is chosen quite arbitrarily in Eqs. ([23]) and (|24|) . To improve the situation, we propose the following procedure to 
determine the S4 dependence of the differential cross sections based on a comparison of the soft-gluon predictions with 
the exact NLO results. First, we define the on-shell Born-level distributions in the LO kinematics, i.e. at S4 = 0, as 



originating from collinear gluon emission. Our numerical analysis shows that the new quantities B^°[ n (s', t±, u\) 
improve essentially the quality of the soft-gluon approximation for both &t and &l- More details can be found in 
Ref. [30(. In our further studies, we use the improved Born-level distributions, B^°[ n (s', t\, u\), instead of old ones 
given by Eqs. (23| an d (|24|) . 

Note that the redefinition of the usual Born-level distributions used in the present paper does not affect any previous 
predictions of the standard resummation approach. The only purpose of our redefinition is to extend the region of 
applicability of the soft-gluon approximation to higher values of Q 2 . 

Our results for the x distributions of the Callan-Gross ratio R(x,Q 2 ) = Fl/Ft in charm leptoproduction are 
presented at several values of £ in the left panel of Fig. [5J For comparison, the K factors K^>(x, Q 2 ) = F^ LO /F^° 
and K^ 2 \x, Q 2 ) = F^ NLO / F^ LO for the transverse structure function are shown at the same values of £ in the right 
panel of Fig. [5] One can see that the sizeable soft-gluon corrections to the production cross sections affect the Born 
predictions for R(x, Q 2 ) both at NLO and NNLO very little, by a few percent only. 

Let us briefly discuss the origin of the perturbative stability of the Callan-Gross ratio. Note that the mere spin- 
independent structure of the Sudakov logarithms can not explain our results, since perturbative stability does not 

(1,0) (0,0) 

take place at the parton level. In fact, the ratios c f 10) (z, Q 2 ) and ^oj(z, Q 2 ) differ essentially from each other, even 
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near the threshold. This is due to the fact that, according to Eq. (j2"2"|) . the soft-gluon corrections are determined 
by convolutions of the Born cross sections with the Sudakov logarithms, which, apart from factorized <5(s4) terms, 
contain also nonfactorizable ones, see Eq. (j2l)|) . For instance, values of z ~ 10 _1 allow s^/m 2 ~ 1 at Q 2 ~ m 2 , which 
leads to significant nonfactorizable corrections. In other words, collinear bremsstrahlung carries away a large part 
of the initial energy. Since the longitudinal and transverse Born-level partonic cross sections have different energy 
behaviors, the so-called soft-gluon radiation has different impacts on these quantities. 

Our analysis shows that two more factors are responsible for perturbative stability of the hadron-level ratio R(x, Q 2 ). 
First, for relatively low virtuality Q 2 ~ m 2 , both &t(z, Q 2 ) and &l{z,Q 2 ) take their maximum values practically at 
the same values of z not far from the threshold. Second, at x ~ 10 _2 -10 _1 , the gluon distribution function supports 
just the threshold region contribution. According to the saddle point arguments, both these factors together lead 
to an approximate factorization of the Sudakov logarithms at the hadron level and essential cancellation of their 
contributions in the ratio R(x, Q 2 ) = Fl/Ft- 

Note also that the situation with perturbative stability of the Callan-Gross ratio is very similar to the corresponding 
one that takes place for the azimuthal asymmetry in heavy-quark photo- and leptoproduction. In detail, the soft-gluon 
corrections to the azimuthal asymmetry were considered in Refs. [lol H3|. 



IV. ANALYTIC LO PREDICTIONS AT LOW x 



Since the radiative corrections to the Callan-Gross ratio in heavy-flavor leptoproduction are small, it makes sense to 
investigate in more detail the corresponding LO predictions. In this Section, we derive compact \ow-x approximation 
formulae for the ratio R2(x,Q 2 ) — IxF^jF^ at LO, which greatly simplify the extraction of the structure function 
Fi(x,Q 2 ) from measurements of the reduced cross section, a(x, Q 2 ), defined by Eqs. J5]) and ©. For this purpose, 
we convolute the LO partonic cross sections given by Eqs. © and (fTUf with the low- a; asymptotics of the gluon PDF: 



1 



(30) 



The value of S in Eq. (|30[) is a matter of discussion. The simplest choice, 5 = 0, leads to a non-singular behavior of the 
structure functions for x — > 0. Another extreme value, 8 = 1/2, historically originates from the BFKL resummation 
of the leading powers of ln(l/a;) In reality, 8 is a function of Q 2 (for an experimental review, see Ref. [3~D]). 

Theoretically, the Q 2 dependence of 8 is calculated using the DGLAP evolution equations [24| . 

First, we calculate the LO hadron-level cross sections for both extreme cases, (5 = and 1/2. Our predictions for 
the quantity R2(x, Q 2 ) in the limit of x — ► have the following form: 



R 



(1/2) 



(Q 2 



2 1 + 6A-4A(1 + 3A)J(A) 
1 + 4A 1 + 2(1- A) J(A) ' 

8 [3 + 4A (13 + 32A)] £(1/(1 + 4A)) - 4A (9 + 32A) K(l/(1 + 4A)) 
1 + 4A (-37 + 72A) £(1/(1 + 4A)) + 2 (23 - 36A) K(l/(l + 4A)) : 



where A is defined in Eq. (fT3 



J(A) 



1 



In 



V1 + 4A+1 



VT + 4A V1 + 4A- 1 

and the functions K(y) and E{y) are the complete elliptic integrals of the first and second kinds defined as 



(31) 
(32) 

(33) 



1 

KW =!7m 



dt 



V(i-* 2 )(i -^ 2 ) 



E(y) = Jdi 



i-yt 2 
l-t 2 



(34) 



The result in Eq. (|3"Tj) was previously found in Ref. [22j], where an approximation to its NLO counterpart was also 
presented. 

The left panel of Fig. [5] shows the ratios R^iQ 2 ) and R^ 2 \Q 2 ) as functions of £. One can see that the difference 
between these quantities varies slowly from 20% at low Q 2 to 10% at high Q 2 . For comparison, also the LO results 
for R2(x,Q 2 ) calculated at several values of x using the CTEQ5L gluon PDF [26| are shown. We observe that, for 
x — > 0, the CTEQ5L predictions converge to the function R% \Q 2 ) practically in the entire region of Q 2 . We have 
verified that the same situation takes also place for all other LO and NLO CTEQ PDF versions [HI, H(| . 
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i;=Q 2 /m 2 ^=Q 2 /m 2 



FIG. 6: LO low- a; predictions for the ratio B,2(x,Q 2 ) — 2xFl/F2 in charm leptoproduction. Left panel: Asymptotic ratios 
R ( 2 0) (Q 2 ) (gray points) and R 2 1/2) (Q 2 ) (black points), as well as CTEQ5L predictions for R 2 {x,Q 2 ) at x = 10~ 2 , 10" 3 and 
10~ 4 . Right panel: Asymptotic ratio R 2 S) (Q 2 ) at S = 0, 0.2, 0.3, 0.4 and 0.5. 



Next, we derive an analytic low- a: formula for the ratio Jig (a;, Q 2 ) with arbitrary values of 5, in terms of the Gauss 



hypergeometric function. Our result has the following form: 



i^(Q 2 ) = 4 



$ ( 1 + 6 > im) (! + 4A ) $ ( 2 + 5 ' TT4A 



1 



5(l-8 2 ) 



(2+<5)(3+5) 

where the function 4>(r, z) is defined as 



(^TTil)-( 1 + 4A ) ( 4 -<5-3^)*( 1 + ^TT4l)' 



(35) 



r(l/2)r(l + r) /l 3 

$(r ' 2) = TT7 r (3 /2 + r) ^^' 1 + r; 2 +r;z 

The hypergeometric function 2 J 7 ! (a, c; z) has the following series expansion: 

r (c) ^ r (a + n) T (6 + n) z n 



2 Fi (a,b;c;z) = 



T(a)T(b) 



E 



T (c + n) 



(36) 



(37) 



In the right panel of Fig. the S dependence of the asymptotic ratio R^iQ 2 ) is investigated. One can see that the 



ratio R { 2 S \q 2 



R { °- 5) (Q 2 ) and R 2 V "'(Q 2 ) varies slowly from 



rapidly converges to the function Jig (Q 2 ) for S > 0.2 



,(0.3) 



In particular, the relative difference between 
i at low Q 2 to 2% at high Q 2 . 
As mentioned above, the Q 2 dependence of the parameter S is determined with the help of the DGLAP evolution. 
However, our analysis shows that hadron-level predictions for i?2(£ — > 0,Q 2 ) depend weakly on S practically in the 
entire region of Q 2 for 0.2 < S < 0.9. For this reason, our simple formula (f3"2")) with 6=1/2 (i.e., without any 
evolution) describes with good accuracy the low- a: CTEQ results for R 2 (x,Q 2 ). We conclude that the hadron-level 
predictions for R 2 (x — * 0, Q 2 ) are stable not only under the NLO corrections to the partonic cross sections, but also 
under the DGLAP evolution of the gluon PDF. 

Finally, we use the analytic expressions (f3"Tj) . (|3l?|) and ([33)1 for the extraction of the structure functions F 2 {x : Q 2 ) 
and F 2 (x, Q 2 ) from the HERA measurements of the reduced cross sections a c (x, Q 2 ) and <r b (x, Q 2 ), respectively. The 
results of our analysis of the HERA data on the charm and bottom electroproduction are collected in Tables H] and 
HTl respectively. In our calculations, the values m c = 1.3 GeV and mi, = 4.3 GeV for the charm and bottom quark 
masses are used. The LO predictions, F 2 (LO), for the cases of S — 0.5, 0.3 and are presented and compared with 
the NLO values, ^(NLO), obtained in the HI analysis [13,[2l[. One can see that all the considered LO predictions 
agree with the NLO results with an accuracy better than 1%. This is because the contributions of the longitudinal 
structure functions, Fl(x,Q 2 ) and F^(x,Q 2 ), to the reduced cross sections, a c (x,Q 2 ) and a b (x, Q 2 ), are small, less 
than 5%, in the kinematic range of the HERA HI experiment. 
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TABLE I: Values of F£(x,Q 2 ) extracted from the HERA measurements of a c (x,Q 2 ) at low [H[ and high [lj Q 2 (in GeV 2 ) 
for various values of x (in units of 10~ 3 ). The NLO HI results [2C| GEO] are compared with the LO predictions corresponding to 
the cases of 5 — 0.5, 0.3 and 0. 



Q 2 


X 


y 


a c 


Error 










(GeV 2 ) (xl0~ 3 ) 






y/o) 


HI 


5 = 0.5 


S = 0.3 


s — o 


12 


0.197 


0.600 


0.412 


18 


0.435 ± 0.078 


0.435 ± 0.078 


0.434 ± 0.078 


0.431 ± 0.077 


12 


0.800 


0.148 


0.185 


13 


0.186 ± 0.024 


f\ -t en K i r\ r\ c\ a 

0.185 ± 0.024 


0.185 ± 0.024 


0.185 ± 0.024 


25 


0.500 


0.492 


0.318 


13 


0.331 ± 0.043 


0.331 ± 0.043 


0.330 ± 0.043 


0.328 ± 0.043 


25 


2.000 


0.123 


0.212 


10 


A Ai A i a A0 1 
U.212 ± 0.021 


a a i o i a a a i 
0.212 ± 0.021 


nmn 1 Anoi 
0.212 ± 0.021 


AOIO 1 AAA1 

0.212 ± 0.021 


60 


2.000 


0.295 


0.364 


10 


0.369 ± 0.040 


0.369 ± 0.040 


0.368 ± 0.040 


0.368 ± 0.040 


60 


5.000 


0.118 


0.200 


12 


0.201 ± 0.024 


0.200 ± 0.024 


0.200 ± 0.024 


0.200 ± 0.024 


200 


0.500 


0.394 


0.197 


23 


0.202 ± 0.046 


0.202 ± 0.046 


0.202 ± 0.046 


0.201 ± 0.046 


200 


1.300 


0.151 


0.130 


24 


0.131 ±0.032 


0.130 ±0.031 


0.130 ±0.031 


0.130 ±0.031 


650 


1.300 


0.492 


0.206 


27 


0.213 ±0.057 


0.213 ±0.057 


0.213 ±0.057 


0.212 ±0.057 


650 


3.200 


0.200 


0.091 


31 


0.092 ± 0.028 


0.091 ±0.028 


0.091 ±0.028 


0.091 ±0.028 



TABLE II: Values of F%(x,Q 2 ) extracted from the HERA measurements of a b (x,Q 2 ) at low [HJ and high [U Q 2 (in GeV 2 ) 
for various values of x (in units of 10~ 3 ). The NLO HI results 0, Hl[ are compared with the LO predictions corresponding to 
the cases of 5 — 0.5, 0.3 and 0. 



Q 2 


X 


y 


d b 


Error 


Fj'(NLO) 


F|(LO) 


F 2 c (LO) 


F|(LO) 


(GeV 2 ) (xHT 3 ) 






(%) 


HI 


S = 0.5 


S = 0.3 


5 = 


12 


0.197 


0.600 


0.0045 


60 


0.0045 ± 0.0027 


0.0046 ± 0.0027 


0.0046 ± 0.0027 


0.0046 ± 0.0027 


12 


0.800 


0.148 


0.0048 


45 


0.0048 ± 0.0022 


0.0048 ± 0.0022 


0.0048 ± 0.0022 


0.0048 ± 0.0022 


25 


0.500 


0.492 


0.0122 


31 


0.0123 ±0.0038 


0.0124 ± 0.0038 


0.0124 ±0.0038 


0.0123 ±0.0038 


25 


2.000 


0.123 


0.0061 


39 


0.0061 ± 0.0024 


0.0061 ± 0.0024 


0.0061 ± 0.0024 


0.0061 ± 0.0024 


60 


2.000 


0.295 


0.0189 


29 


0.0190 ± 0.0055 


0.0190 ± 0.0055 


0.0190 ±0.0055 


0.0190 ±0.0055 


60 


5.000 


0.118 


0.0130 


36 


0.0130 ± 0.0047 


0.0130 ± 0.0047 


0.0130 ± 0.0047 


0.0130 ± 0.0047 


200 


0.500 


0.394 


0.0393 


31 


0.0413 ±0.0128 


0.0402 ± 0.0125 


0.0401 ± 0.0125 


0.0400 ± 0.0124 


200 


1.300 


0.151 


0.0212 


38 


0.0214 ±0.0081 


0.0213 ± 0.0081 


0.0213 ±0.0081 


0.0212 ±0.0081 


650 


1.300 


0.492 


0.0230 


51 


0.0243 ± 0.0124 


0.0240 ±0.0122 


0.0239 ± 0.0122 


0.0238 ± 0.0121 


650 


3.200 


0.200 


0.0124 


44 


0.0125 ±0.0055 


0.0125 ±0.0055 


0.0125 ±0.0055 


0.0125 ±0.0055 



V. CONCLUSION 

We conclude by summarizing our main observations. In the present paper, we studied the radiative corrections 
to the Callan-Gross ratio R(x, Q 2 ) in heavy-quark leptoproduction. We considered the exact NLO results at low 
and moderate Q 2 < in 2 , asymptotic NLO predictions at high Q 2 ^S> m 2 , and both NLO and NNLO soft-gluon (or 
threshold) corrections at large Bjorken x. It turned out that large (especially, at non-small x) radiative corrections 
to the structure functions Ft(x,Q 2 ) and Fl(x,Q 2 ) cancel each other in their ratio R(x,Q 2 ) — Fl/Ft with good 
accuracy. As a result, the NLO contributions to the ratio R(x, Q 2 ) are less than 10% in a wide region of the variables 
x and Q 2 . Our analysis also shows that the NLO predictions for R(x,Q 2 ) are sufficiently insensitive (to within 
ten percent) to standard uncertainties in the QCD input parameters. We conclude that, unlike the production cross 
sections, the Callan-Gross ratio in heavy-quark leptoproduction is quantitatively well defined in pQCD. Measurements 
of the quantity R(x, Q 2 ) in charm and bottom leptoproduction would provide a good test of the conventional parton 
model based on pQCD. 

Concerning the experimental aspects, we propose to exploit the observed perturbative stability of the Callan-Gross 
ratio in the extraction of the structure functions F£(x, Q 2 ) and F$(x, Q 2 ) from the corresponding reduced cross 
sections. For this purpose, we provided compact LO hadron-level formulae for the ratio R2{x,Q 2 ) = 2xFl/F 2 — 
R/(l + R) in the limit x — > 0. We demonstrated that these analytic expressions simplify the extraction of i^x, Q 2 ) 
without affecting the accuracy of the result in practice. In particular, our LO formula for R2(x, Q 2 ) with 8 = 1/2 
usefully reproduces the results for F^x^Q 2 ) and F!l(x,Q 2 ) obtained by the HI Collaboration [2fJ,[2l| with the help 
of the more cumbersome NLO evaluation of Fl(x, Q 2 ). 
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In this paper, we investigated the contribution to R{x,Q 2 ) from the dominant mechanism, photon-gluon fusion, 
within the fixed-flavor-number scheme. To take into account the photon-heavy-quark scattering component, one 
should adopt the variable-flavor-number scheme, which allows one to rcsum potentially large mass logarithms of the 
type a s In (Q 2 /m 2 ) , whose contribution dominates at Q 2 ^S> m 2 . Some recent developments concerning this scheme 
may be found in Ref. [32 ]. The variable-flavor-number-scheme predictions for the Callan-Gross ratio as well as the 
possibility to discriminate experimentally between photon-gluon fusion and quark-scattering contributions to R(x, Q 2 ) 
will be considered in a forthcoming publication. 
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